Order of Two-dimensional Isotropic 
Dipolar Antiferromagnets 

C. PICH and F. SCHWABL 

Institut fur Theoretische Physik 
Physik- Department der Technischen Universitdt Miinchen 
D-8046 Garching, Federal Republic of Germany 

The question of the existence of order in two-dimensional isotropic dipolar Heisen- 
berg antiferromagnets is studied. It is shown that the dipolar interaction leads to 
a gap in the spin- wave energy and a nonvanishing order parameter. The resulting 
finite Neel-temperature is calculated for a square lattice by means of linear spin- wave 
theory. 

PACS numbers: 75.10 J, 75.30 D, 75.30 G, 75.50 E 



Typeset Using REVTEX 



1 



The question of order in low-dimensional systems has attracted the interest of theoretical 
and experimental physics for a long time. It has been pointed out already by Bloch and 
has been proven exactly @ that long-range order is absent in isotropic two-dimensional 
Heisenberg ferromagnets with short-range interaction. The same is true for one- and two- 
dimensional He 4 and antiferromagnets. 

In real systems one unavoidably has a dipolar interaction in addition to the short-range 
exchange interaction, which breaks the rotational symmetry. It has been shown by Maleev 
f| that the q 2 dispersion law of the isotropic ferromagnet is modified such that a finite order 
parameter exists in two dimensions. The finite temperature behaviour and in particular the 
transition temperature have been calculated by Pokrovsky and Feigelman ||. 

At first glance one might believe that the dipolar interaction is unimportant in antiferro- 
magnets due to cancellations because of the alternating order and it may come as a surprise 
that this expectation is incorrect. A first hint that the dipolar interaction can influence 
antiferromagnetic behaviour comes from the critical region. There the nonlinear coupling of 
fluctuations of the staggered magnetization and the magnetization, which is no more con- 
served, leads to a change of the critical dynamic exponent and of the scaling functions [0. 
In the low temperature phase antiferromagnetic spin-waves involve the coupled precessional 
motion of magnetization and staggered magnetization. Since the conservation law for the 
magnetization is broken by the dipolar forces also the magnon frequency becomes finite at 
wave-vector q = 0. Thus we will show below that (i) two-dimensional Heisenberg antiferro- 
magnetic order exists on simple square lattices due to the dipolar interaction, with the spin 
orientation perpendicular to the plane, (ii) the magnon frequency has a gap, the magnitude 
of which depends on the dipolar interaction and the exchange interaction, (iii) the critical 
field for which the spins rearrange is finite, (iv) there is a finite Neel-temperature which is 
evaluated. 

The Hamiltonian of a dipolar antiferromagnet reads 

H = - £ £ (>"W + A$) S?S? r - g^ B H £ Sf (1) 
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with spins Si at lattice sites x/. The first term in brackets is the exchange interaction Ju> 
and the second the dipole-dipole interaction with 

A*? = -V*) 2 ( , 5ap „ - m^p^Lz^A . (2) 

2 \ — x Z / 1 3 \xi-x v \ b J 

Although we are mainly interested in these two terms we have also included a homogeneous 
external field H along the z-axis (g denotes the Lande factor and hb the Bohr magneton). In 
a two-dimensional system with additional dipole-dipole interaction the rotational symmetry 
is broken; thus the theorem of Hohenberg-Mermin- Wagner 0|| does not apply. 

We consider a square lattice in the xy plane with lattice constant a and the spins orien- 
tated alternatingly along the z axis. By means of the Holstein-Primakoff transformation || 
the Hamilton operator can be expressed in terms of the Bose operators (aj, a{) (neglecting 
terms higher than bilinear in Eqs. |^ and |j) 

S? = ^{o l + aS) t Sf = Ti^(ai - a?) , Sf = ±(S - o^o,), (3) 

where the upper (lower) sign is for the first (second) sublattice. This transformation and a 
Fourier transformation yields 

H = T,<r{A$ aiciq- + ^B^ (a ? a_ ? + ajat^) + 



C 9 aa a_ ? -_ g - + <7* a J-« T _-_-. + D 9 aJ-« 9 - +(fo + D* a} +$) a ? } (4) 

with the coefficients 

Aq = S(2J$ — Jq — Jq+q ) + S(2Ag* — A~ x — Af^ Q ) (5a) 

Bf = S(J €+€o - Jq) + S(A^ o - Af ) (5b) 



Q = iSAf (5c) 
1 

T 



D ? =iSAf + -gpL B H . (5d) 



In this description (Eq. ||) the primitive cell is the chemical 0, which is half the magnetic. 
The wave vector go = ^ (1, 1) represents the antiferromagnetic, staggered modulation via 
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e iqoxi _ rj^g are p our j er transform of the dipole tensor (Eq. |2|) and can be calculated 



by the method of Ewald summation [10 



The Hamiltonian (Eq. |j) is diagonalized by a generalized Bogoliubov transformation 
with two kinds of creation and annihilation operators c^, , c|, c| 

#=£(o)+£]T44 + 4 (6) 
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=£44+ 4 c -? + 4+?o4+ 9 -b + ^A-^ ( 7 ) 

i=l 

4>4 t l = ( W^> 4,ci] = [cj/, c /] = o, 

with wave-vectors restricted to the magnetic Brillouin zone. Here E(0) is the ground-state 
energy. The spin-wave energies then assume the form 

4 a = ±(n! ± n a ) (8) 

with 
and 

^2 = ~ B\~ Aq+qo + -^l+gi,) 2 + 16 [C^, ( A^ Q — -Bg+g ) — C$ (Aj — B^)\ 

X [Cq (Acf+ql + B$ + f ) — Ctf + tf Q (Af+ B^j] 

+4( W tf ) 2 ((^+ ^ +(fo ) 2 - (Bj-B^) 2 ). 

Let us now discuss Eq. (8) in the case of primary interest namely vanishing external field 
(H = 0). The dipolar interaction has two effects: first, in contrast to the isotropic case, 
the excitation spectrum is no more degenerate, i.e. two different branches appear. Second 
it produces an energy gap for q — > 



E = 2S^A% - A^(J ql] - Jo) - « - A%) (9) 

with 

A f = \( A f + A F)- 

In Fig. I] the dispersion relation is shown for three values for the ratio of dipolar and 
exchange energy k = ^fjr^s with isotropic nearest-neighbor exchange interaction (J < 0). 
The two branches can be resolved only for large values of k. For more realistic ratios (10~ 3 ) 
the two magnon branches practically coincide, but a significant deviation from the pure 
exchange case still remains in the immediate vicinity of the zone center. The argument 
of the first square root in Eq. @ for the gap equals the difference of dipolar energy for 
out-of- and in-plane staggered orientation and is positive. Thus stability of the ground state 
requires 

j a - Jo > K p - A£ > o. 

For dipole interaction alone the system would order with the magnetization in the plane. 
Thus the dipolar energy difference between in-plane and out-of plane orientation must be 
exceeded by the exchange energy in order to favour the assumed configuration. In particular 
the gap is proportional to the square root of the difference of the static energy between the 
configurations of in-plane and out-of-plane magnetization. In a three-dimensional simple cu- 
bic lattice the first root in Eq. @ vanishes because of the symmetry, but in two-dimensional 
systems there is a finite gap for perpendicular antiferromagnetic order. We note that for 
sufficiently large exchange energy the gap is the geometric mean of dipole and exchange 
energy, which in turn implies that the gap is much larger than the dipolar energy for k < 1. 

Let us add some comments on the interplay of exchange and dipolar interaction. The 
former imposes the antiferromagnetic order while the latter leads to the orientation perpen- 
dicular to the plane and prevents thermal fluctuations from its destruction. To exhibit more 
clearly the physical origin of the energy gap and its principal dependence on dipolar and 
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exchange interaction we exhibit the equations of motion for the spin components. Approxi- 
mating the longitudinal part (z— component) by Sf ~ Se iq ° Xl and specializing to q — 0, the 
equations for the transverse components become 

SS= (Ao-B )Sl (10a) 

5J = - (A + B )S* (10b) 

with an analogous set for the Sq component. The coefficient on the right hand side of 
Eq. (|10a|) assumes a finite value in contrast to pure exchange antiferromagnets where 5| is 
conserved. Thus the coupled motion of Sq and S^ o leads precisely to the finite energy gap 



E (Eq. ©) [|TT 



From the spin- wave energy (Eq. @) we can calculate the critical field for which the 
antiferromagnetic Neel state gets destabilized by a magnetic field. It is given by the field 
Hq for which the energy (q = 0) vanishes 

H*=— E . (11) 

Hence the critical field is proportional to the energy gap. In 2D the anisotropic dipolar 
interaction stabilizes the antiferromagnetic configuration in an external field up to the above 
value. 

Now we turn to the evaluation of T/v the transition temperature for vanishing external 
field, i.e. the temperature at which the staggered magnetization vanishes. We use linear 
spin-wave theory, i.e. interactions between magnons and temperature renormalization of 
the magnon energy are neglected. This approximation is justified at low temperature and 
should lead to an order of magnitude estimate of the main dependence on exchange and 
dipolar interaction. The staggered magnetization then reads 

N(T) = gf i B (NS -£< a\ ai >) = m - N - N th (T). (12) 
i 

This sum is calculated by means of the transformation (Eq. [?[) and requires the evaluation 
of the coefficients u q , v q , s q , t q , which are complicated functions of the the coefficients in Eqs. 
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(5a-5d). The number of thermally excited magnons N t h{T) can be expressed in terms of the 
mean number of excitations 

4=<4 t 4>=( e ^ fcflT - 1 )" 1 - ( 13 ) 

We now consider isotropic nearest-neighbor exchange (J < 0) and focuse only on the limit 
of small wave- vectors and small dipolar energies (k <C 1). The deviation of the ground-state 
magnetization due to thermal excitations then takes the form 

N*(T) = (14) 

and the deviation originating from the zero-point oscillations 

No = ^E^-l) (15) 



with D = 8S\J\ and Eg « J D 2 a 2 q 2 /2 + Eq. The zero point deviation N Q has a finite value 
ill for pure exchange interaction already, and is effected only neglegibly by the dipole inter- 



action. The dipole interaction favours the antiferromagnetic order and leads to a reduction 
of N of order k. 

Now we turn to N t h{T), which is divergent for pure exchange antiferromagnets implying 
the absence of antiferromagnetic order in this case. The small wave-vector approximation 
in Eqs. (13-15) is not accurate for temperatures near the phase transition but it is sufficient 
for our crude estimate. From Fig. H it becomes clear that the regime of small wave-vectors 
is essential for the calculation of the sum. The existence of a nonvanishing gap makes the 
sum convergent and allows a phase-transition (Eq. |12|) at a finite temperature TV 

N th (T N ) = m - N = b'. (16) 

After replacing the sum in Eq. (14) by an integral and the Brillouin zone by a circle of the 
same area, the evaluation leads to the following implicit equation (the upper bound of the 
integral is set to infinity) 



for T/v with b = 2 ^n ^'' ^ n ^ ne °f vanishing gap (Eq — > 0) we recover again the 

impossibility of a phase-transition. From Eq. (|i~7D an asymptotic solution for small dipole 
energies can be derived 



r " = -^W~i4' (18) 



We now compare our results with experiments on K^MnF^ for which the spin- wave 
dispersion has been measured K 2 MnF^ is a quasi two-dimensional antiferromagnet 

with the spin orientation perpendicular to the ab-plane, a transition temperature T/v = 42 K 
and an exchange energy | J\\ = 8.7 K (S = 5/2). An energy gap of Eq = 7.5 K is observed. 
Evaluation of the energy gap via Eq. (9) yields Eq = 7.6 K, which is in remarkably agreement 
with the experimental result. Solving Eq. (17) with the zero point deviation iVo = 0.2Ag/i£ 
|12]| our estimate for the transition temperature becomes T/v = 112 K which is too large 
by a factor of three. It will be lowered if instead of the small wave vector expansion the 
correct dispersion relation (Eq. 8) shown in Fig. 1 is used in the evaluation of Nth (Eq. 14). 
Furthermore we have neglected the interaction of magnons, which will be important at higher 



temperatures and will lower T/v- This could be treated by more elaborate theories, e.g. |T4 
but goes beyond the scope of this paper. 

In summary we have shown that two-dimensional antiferromagnetic order is possible due 
to the dipolar interaction. 
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FIGURES 

FIG. 1. The spin- wave dispersion relation (Eq. ||) of pure exchange antiferromagnets with 
nearest-neighbour interaction (solid line) and with additional dipolar interaction (S = 1/2), for 
the ratios of dipolar energy to exchange energy k = ^rjra along the 0] direction: k = 0.1 

(dashed), k = 0.01 (dot-dashed) and k = 0.001 (dot-dashed-dashed). The splitting of the two 
magnon branches is visible only for k = 0.1. 
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